Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

RNA sequencing unravels novel L cell constituents and mechanisms of GLP-1 secretion in human gastric bypass-operated intestine.

AIMS/HYPOTHESIS: Roux-en-Y gastric bypass surgery (RYGB) frequently results in remission of type 2 diabetes as well as exaggerated secretion of glucagon-like peptide-1 (GLP-1). Here, we assessed RYGB-induced transcriptomic alterations in the small intestine and investigated how they were related to the regulation of GLP-1 production and secretion in vitro and in vivo. METHODS: Human jejunal samples taken perisurgically and 1 year post RYGB (n=13) were analysed by RNA-seq. Guided by bioinformatics analysis we targeted four genes involved in cholesterol biosynthesis, which we confirmed to be expressed in human L cells, for potential involvement in GLP-1 regulation using siRNAs in GLUTag and STC-1 cells. Gene expression analyses, GLP-1 secretion measurements, intracellular calcium imaging and RNA-seq were performed in vitro. OGTTs were performed in C57BL/6j and iScd1-/- mice and immunohistochemistry and gene expression analyses were performed ex vivo. RESULTS: Gene Ontology (GO) analysis identified cholesterol biosynthesis as being most affected by RYGB. Silencing or chemical inhibition of stearoyl-CoA desaturase 1 (SCD1), a key enzyme in the synthesis of monounsaturated fatty acids, was found to reduce Gcg expression and secretion of GLP-1 by GLUTag and STC-1 cells. Scd1 knockdown also reduced intracellular Ca2+ signalling and membrane depolarisation. Furthermore, Scd1 mRNA expression was found to be regulated by NEFAs but not glucose. RNA-seq of SCD1 inhibitor-treated GLUTag cells identified altered expression of genes implicated in ATP generation and glycolysis. Finally, gene expression and immunohistochemical analysis of the jejunum of the intestine-specific Scd1 knockout mouse model, iScd1-/-, revealed a twofold higher L cell density and a twofold increase in Gcg mRNA expression. CONCLUSIONS/INTERPRETATION: RYGB caused robust alterations in the jejunal transcriptome, with genes involved in cholesterol biosynthesis being most affected. Our data highlight SCD as an RYGB-regulated L cell constituent that regulates the production and secretion of GLP-1.

GPR162 is a beta cell CART receptor.

Cocaine and amphetamine-regulated transcript (CART) is expressed in pancreatic islet cells and neuronal elements. We have previously established insulinotropic actions of CART in human and rodent islets. The receptor for CART in the pancreatic beta cells is unidentified. We used RNA sequencing of Cartpt knockdown (KD) INS-1 832/13 cells and identified GPR162 as the most Cartpt-regulated receptor. We therefore tested if GPR162 mediates the effects of CART in beta cells. Binding of CART to GPR162 was established using proximity ligation assay, radioactive binding, and co-immunoprecipitation, and KD of Gpr162 mRNA caused reduced binding. Gpr162 KD cells had blunted CARTp-induced exocytosis, and reduced CARTp-induced insulin secretion. Furthermore, we identified a hitherto undescribed GPR162-dependent role of CART as a regulator of cytoskeletal arrangement. Thus, our findings provide mechanistic insight into the effect of CART on insulin secretion and show that GPR162 is the CART receptor in beta cells.

How broad is the selfing syndrome? Insights from convergent evolution of gene expression across species and tissues in the Capsella genus.

The shift from outcrossing to selfing is one of the main evolutionary transitions in plants. It is accompanied by profound effects on reproductive traits, the so-called selfing syndrome. Because the transition to selfing also implies deep genomic and ecological changes, one also expects to observe a genomic selfing syndrome. We took advantage of the three independent transitions from outcrossing to selfing in the Capsella genus to characterize the overall impact of mating system change on RNA expression, in flowers but also in leaves and roots. We quantified the extent of both selfing and genomic syndromes, and tested whether changes in expression corresponded to adaptation to selfing or to relaxed selection on traits that were constrained in outcrossers. Mating system change affected gene expression in all three tissues but more so in flowers than in roots and leaves. Gene expression in selfing species tended to converge in flowers but diverged in the two other tissues. Hence, convergent adaptation to selfing dominates in flowers, whereas genetic drift plays a more important role in leaves and roots. The effect of mating system transition is not limited to reproductive tissues and corresponds to both adaptation to selfing and relaxed selection on previously constrained traits.

Myocardin regulates exon usage in smooth muscle cells through induction of splicing regulatory factors.

Differentiation of smooth muscle cells (SMCs) depends on serum response factor (SRF) and its co-activator myocardin (MYOCD). The role of MYOCD for the SMC program of gene transcription is well established. In contrast, the role of MYOCD in control of SMC-specific alternative exon usage, including exon splicing, has not been explored. In the current work we identified four splicing factors (MBNL1, RBPMS, RBPMS2, and RBFOX2) that correlate with MYOCD across human SMC tissues. Forced expression of MYOCD family members in human coronary artery SMCs in vitro upregulated expression of these splicing factors. For global profiling of transcript diversity, we performed RNA-sequencing after MYOCD transduction. We analyzed alternative transcripts with three different methods. Exon-based analysis identified 1637 features with differential exon usage. For example, usage of 3´ exons in MYLK that encode telokin increased relative to 5´ exons, as did the 17 kDa telokin to 130 kDa MYLK protein ratio. Dedicated event-based analysis identified 239 MYOCD-driven splicing events. Events involving MBNL1, MCAM, and ACTN1 were among the most prominent, and this was confirmed using variant-specific PCR analyses. In support of a role for RBPMS and RBFOX2 in MYOCD-driven splicing we found enrichment of their binding motifs around differentially spliced exons. Moreover, knockdown of either RBPMS or RBFOX2 antagonized splicing events stimulated by MYOCD, including those involving ACTN1, VCL, and MBNL1. Supporting an in vivo role of MYOCD-SRF-driven splicing, we demonstrate altered Rbpms expression and splicing in inducible and SMC-specific Srf knockout mice. We conclude that MYOCD-SRF, in part via RBPMS and RBFOX2, induce a program of differential exon usage and alternative splicing as part of the broader program of SMC differentiation.

High-throughput muscle fiber typing from RNA sequencing data.

BACKGROUND: Skeletal muscle fiber type distribution has implications for human health, muscle function, and performance. This knowledge has been gathered using labor-intensive and costly methodology that limited these studies. Here, we present a method based on muscle tissue RNA sequencing data (totRNAseq) to estimate the distribution of skeletal muscle fiber types from frozen human samples, allowing for a larger number of individuals to be tested. METHODS: By using single-nuclei RNA sequencing (snRNAseq) data as a reference, cluster expression signatures were produced by averaging gene expression of cluster gene markers and then applying these to totRNAseq data and inferring muscle fiber nuclei type via linear matrix decomposition. This estimate was then compared with fiber type distribution measured by ATPase staining or myosin heavy chain protein isoform distribution of 62 muscle samples in two independent cohorts (n = 39 and 22). RESULTS: The correlation between the sequencing-based method and the other two were rATPas = 0.44 [0.13-0.67], [95% CI], and rmyosin = 0.83 [0.61-0.93], with p = 5.70 × 10-3 and 2.00 × 10-6, respectively. The deconvolution inference of fiber type composition was accurate even for very low totRNAseq sequencing depths, i.e., down to an average of ~ 10,000 paired-end reads. CONCLUSIONS: This new method ( https://github.com/OlaHanssonLab/PredictFiberType ) consequently allows for measurement of fiber type distribution of a larger number of samples using totRNAseq in a cost and labor-efficient way. It is now feasible to study the association between fiber type distribution and e.g. health outcomes in large well-powered studies.

The relative role of plasticity and demographic history in Capsella bursa-pastoris: a common garden experiment in Asia and Europe.

The colonization success of a species depends on the interplay between its phenotypic plasticity, adaptive potential and demographic history. Assessing their relative contributions during the different phases of a species range expansion is challenging, and requires large-scale experiments. Here, we investigated the relative contributions of plasticity, performance and demographic history to the worldwide expansion of the shepherd's purse, Capsella bursa-pastoris. We installed two large common gardens of the shepherd's purse, a young, self-fertilizing, allopolyploid weed with a worldwide distribution. One common garden was located in Europe, the other in Asia. We used accessions from three distinct genetic clusters (Middle East, Europe and Asia) that reflect the demographic history of the species. Several life-history traits were measured. To explain the phenotypic variation between and within genetic clusters, we analysed the effects of (i) the genetic clusters, (ii) the phenotypic plasticity and its association to fitness and (iii) the distance in terms of bioclimatic variables between the sampling site of an accession and the common garden, i.e. the environmental distance. Our experiment showed that (i) the performance of C. bursa-pastoris is closely related to its high phenotypic plasticity; (ii) within a common garden, genetic cluster was a main determinant of phenotypic differences; and (iii) at the scale of the experiment, the effect of environmental distance to the common garden could not be distinguished from that of genetic clusters. Phenotypic plasticity and demographic history both play important role at different stages of range expansion. The success of the worldwide expansion of C. bursa-pastoris was undoubtedly influenced by its strong phenotypic plasticity.

Correction: The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels.

[This corrects the article DOI: 10.1371/journal.pgen.1009726.].

The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels.

Selective breeding for desirable traits in strictly controlled populations has generated an extraordinary diversity in canine morphology and behaviour, but has also led to loss of genetic variation and random entrapment of disease alleles. As a consequence, specific diseases are now prevalent in certain breeds, but whether the recent breeding practice led to an overall increase in genetic load remains unclear. Here we generate whole genome sequencing (WGS) data from 20 dogs per breed from eight breeds and document a ~10% rise in the number of derived alleles per genome at evolutionarily conserved sites in the heavily bottlenecked cavalier King Charles spaniel breed (cKCs) relative to in most breeds studied here. Our finding represents the first clear indication of a relative increase in levels of deleterious genetic variation in a specific breed, arguing that recent breeding practices probably were associated with an accumulation of genetic load in dogs. We then use the WGS data to identify candidate risk alleles for the most common cause for veterinary care in cKCs-the heart disease myxomatous mitral valve disease (MMVD). We verify a potential link to MMVD for candidate variants near the heart specific NEBL gene in a dachshund population and show that two of the NEBL candidate variants have regulatory potential in heart-derived cell lines and are associated with reduced NEBL isoform nebulette expression in papillary muscle (but not in mitral valve, nor in left ventricular wall). Alleles linked to reduced nebulette expression may hence predispose cKCs and other breeds to MMVD via loss of papillary muscle integrity.

Regulation of the Muscarinic M3 Receptor by Myocardin-Related Transcription Factors.

Myocardin-related transcription factors (MRTFs: myocardin/MYOCD, MRTF-A/MRTFA, and MRTF-B/MRTFB) are co-factors of serum response factor (SRF) that activate the smooth muscle cell (SMC) gene program and that play roles in cardiovascular development and mechanobiology. Gain and loss of function experiments have defined the SMC gene program under control of MRTFs, yet full understanding of their impact is lacking. In the present study, we tested the hypothesis that the muscarinic M3 receptor (CHRM3) is regulated by MRTFs together with SRF. Forced expression of MYOCD (8d) in human coronary artery (SMC) followed by RNA-sequencing showed increased levels of M2, M3, and M5 receptors (CHRM2: 2-fold, CHRM3: 16-fold, and CHRM5: 2-fold). The effect of MYOCD on M3 was confirmed by RT-qPCR using both coronary artery and urinary bladder SMCs, and correlation analyses using human transcriptomic datasets suggested that M3 may also be regulated by MRTF-B. Head-to-head comparisons of MYOCD, MRTF-A and MRTF-B, argued that while all MRTFs are effective, MRTF-B is the most powerful transactivator of CHRM3, causing a 600-fold increase at 120h. Accordingly, MRTF-B conferred responsiveness to the muscarinic agonist carbachol in Ca2+ imaging experiments. M3 was suppressed on treatment with the MRTF-SRF inhibitor CCG-1423 using SMCs transduced with either MRTF-A or MRTF-B and using intact mouse esophagus in culture (by 92±2%). Moreover, silencing of SRF with a short hairpin reduced CHRM3 (by >60%) in parallel with α-actin (ACTA2). Tamoxifen inducible knockout of Srf in smooth muscle reduced Srf (by 54±4%) and Chrm3 (by 41±6%) in the urinary bladder at 10days, but Srf was much less reduced or unchanged in aorta, ileum, colon, trachea, and esophagus. Longer induction (21d) further accentuated the reduction of Chrm3 in the bladder and ileum, but no change was seen in the aorta. Single cell RNA-sequencing revealed that Mrtfb dominates in ECs, while Myocd dominates in SMCs, raising the possibility that Chrm3 may be driven by Mrtfb-Srf in the endothelium and by Myocd-Srf in SMCs. These findings define a novel transcriptional control mechanism for muscarinic M3 receptors in human cells, and in mice, that could be targeted for therapy.

Fifteen shades of green: The evolution of Bufotes toads revisited.

The radiation of Palearctic green toads (Bufotes) holds great potential to evaluate the role of hybridization in phylogeography at multiple stages along the speciation continuum. With fifteen species representing three ploidy levels, this model system is particularly attractive to examine the causes and consequences of allopolyploidization, a prevalent yet enigmatic pathway towards hybrid speciation. Despite substantial efforts, the evolutionary history of this species complex remains largely blurred by the lack of consistency among the corresponding literature. To get a fresh, comprehensive view on Bufotes phylogeography, here we combined genome-wide multilocus analyses (RAD-seq) with an extensive compilation of mitochondrial, genome size, niche modelling, distribution and phenotypic (bioacoustics, morphometrics, toxin composition) datasets, representing hundreds of populations throughout Eurasia. We provide a fully resolved nuclear phylogeny for Bufotes and highlight exceptional cyto-nuclear discordances characteristic of complete mtDNA replacement (in 20% of species), mitochondrial surfing during post-glacial expansions, and the formation of homoploid hybrid populations. Moreover, we traced the origin of several allopolyploids down to species level, showing that all were exclusively fathered by the West Himalayan B. latastii but mothered by several diploid forms inhabiting Central Asian lowlands, an asymmetry consistent with hypotheses on mate choice and Dobzhansky-Muller incompatibilities. Their intermediate call phenotypes potentially allowed for rapid reproductive isolation, while toxin compositions converged towards the ecologically-closest parent. Across the radiation, we pinpoint a stepwise progression of reproductive isolation through time, with a threshold below which hybridizability is irrespective of divergence (<6My), above which species barely admix and eventually evolve different mating calls (6-10My), or can successfully cross-breed through allopolyploidization (>15My). Finally, we clarified the taxonomy of Bufotes (including genetic analyses of type series) and formally described two new species, B. cypriensis sp. nov. (endemic to Cyprus) and B. perrini sp. nov. (endemic to Central Asia). Embracing the genomic age, our framework marks the advent of a new exciting era for evolutionary research in these iconic amphibians.

Towards the new normal: Transcriptomic convergence and genomic legacy of the two subgenomes of an allopolyploid weed (Capsella bursa-pastoris).

Allopolyploidy has played a major role in plant evolution but its impact on genome diversity and expression patterns remains to be understood. Some studies found important genomic and transcriptomic changes in allopolyploids, whereas others detected a strong parental legacy and more subtle changes. The allotetraploid C. bursa-pastoris originated around 100,000 years ago and one could expect the genetic polymorphism of the two subgenomes to follow similar trajectories and their transcriptomes to start functioning together. To test this hypothesis, we sequenced the genomes and the transcriptomes (three tissues) of allotetraploid C. bursa-pastoris and its parental species, the outcrossing C. grandiflora and the self-fertilizing C. orientalis. Comparison of the divergence in expression between subgenomes, on the one hand, and divergence in expression between the parental species, on the other hand, indicated a strong parental legacy with a majority of genes exhibiting a conserved pattern and cis-regulation. However, a large proportion of the genes that were differentially expressed between the two subgenomes, were also under trans-regulation reflecting the establishment of a new regulatory pattern. Parental dominance varied among tissues: expression in flowers was closer to that of C. orientalis and expression in root and leaf to that of C. grandiflora. Since deleterious mutations accumulated preferentially on the C. orientalis subgenome, the bias in expression towards C. orientalis observed in flowers indicates that expression changes could be adaptive and related to the selfing syndrome, while biases in the roots and leaves towards the C. grandiflora subgenome may be reflective of the differential genetic load.

Parental legacy, demography, and admixture influenced the evolution of the two subgenomes of the tetraploid Capsella bursa-pastoris (Brassicaceae).

Allopolyploidy is generally perceived as a major source of evolutionary novelties and as an instantaneous way to create isolation barriers. However, we do not have a clear understanding of how two subgenomes evolve and interact once they have fused in an allopolyploid species nor how isolated they are from their relatives. Here, we address these questions by analyzing genomic and transcriptomic data of allotetraploid Capsella bursa-pastoris in three differentiated populations, Asia, Europe, and the Middle East. We phased the two subgenomes, one descended from the outcrossing and highly diverse Capsella grandiflora (CbpCg) and the other one from the selfing and genetically depauperate Capsella orientalis (CbpCo). For each subgenome, we assessed its relationship with the diploid relatives, temporal changes of effective population size (Ne), signatures of positive and negative selection, and gene expression patterns. In all three regions, Ne of the two subgenomes decreased gradually over time and the CbpCo subgenome accumulated more deleterious changes than CbpCg. There were signs of widespread admixture between C. bursa-pastoris and its diploid relatives. The two subgenomes were impacted differentially depending on geographic region suggesting either strong interploidy gene flow or multiple origins of C. bursa-pastoris. Selective sweeps were more common on the CbpCg subgenome in Europe and the Middle East, and on the CbpCo subgenome in Asia. In contrast, differences in expression were limited with the CbpCg subgenome slightly more expressed than CbpCo in Europe and the Middle-East. In summary, after more than 100,000 generations of co-existence, the two subgenomes of C. bursa-pastoris still retained a strong signature of parental legacy but their evolutionary trajectory strongly varied across geographic regions.

The influence of population structure on gene expression and flowering time variation in the ubiquitous weed Capsella bursa-pastoris (Brassicaceae).

Population structure is a potential problem when testing for adaptive phenotypic differences among populations. The observed phenotypic differences among populations can simply be due to genetic drift, and if the genetic distance between them is not considered, the differentiation may be falsely interpreted as adaptive. Conversely, adaptive and demographic processes might have been tightly associated and correcting for the population structure may lead to false negatives. Here, we evaluated this problem in the cosmopolitan weed Capsella bursa-pastoris. We used RNA-Seq to analyse gene expression differences among 24 accessions, which belonged to a much larger group that had been previously characterized for flowering time and circadian rhythm and were genotyped using genotyping-by-sequencing (GBS) technique. We found that clustering of accessions for gene expression retrieved the same three clusters that were obtained with GBS data previously, namely Europe, the Middle East and Asia. Moreover, the three groups were also differentiated for both flowering time and circadian rhythm variation. Correction for population genetic structure when analysing differential gene expression analysis removed all differences among the three groups. This may suggest that most differences are neutral and simply reflect population history. However, geographical variation in flowering time and circadian rhythm indicated that the distribution of adaptive traits might be confounded by population structure. To bypass this confounding effect, we compared gene expression differentiation between flowering ecotypes within the genetic groups. Among the differentially expressed genes, FLOWERING LOCUS C was the strongest candidate for local adaptation in regulation of flowering time.